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ABSTRACT 

We present three XMM-Newton observations of the ultra-luminous compact X-ray source Holmberg 
II X-l in its historical brightest and faintest states. The source was in its brightest state in April 2002 
with an isotropic X-ray luminosity of ~ 2 x 10 40 erg s _1 but changed to a peculiar low/soft state in 
September 2002 in which the X-ray flux dropped by a factor of ~ 4 and the spectrum softened. In all 
cases, a soft excess component, which can be described by a simple or multicolor disk blackbody (MCD; 
kT ~ 120 — 170 eV), is statistically required in addition to a power-law continuum (T ~ 2.4 — 2.9). Both 
the spectral components became weaker and softer in the low/soft state, however, the dramatic variability 
is seen in the power-law component. This spectral transition is opposite to the 'canonical' high/soft - 
low/hard transitions seen in many Galactic black hole binaries. There is possible contribution from 
an optically thin thermal plasma. When this component is taken into account, the spectral transition 
appears to be normal - a drop of the power-law flux and slightly softer blackbody component in the low 
state. 

Subject headings: accretion disks — stars: individual (Holmberg II X-l) — X-rays: stars 



1. INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are the compact, 
off-nuclear galactic X-ray sources with isotropic lumi- 
nosity exceeding the Eddington limit for a stellar mass 
( ^ 2OA/0) black hole (BH). Isotropic emission from in- 
termediate mass (~ 20 — 10 5 M©) BHs and anisotropic 
emission from stellar mass BHs are both attractive inter- 
pretation for the nature of ULXs, however, there is no defi- 
nite proof for either of the two scenarios (Miller & Colbert 
2003). Several observations suggest that ULXs may be 
similar to the X-ray binaries. Discovery of orbital modu- 
lations from several ULXs (Bauer et al. 2001; Sugiho et al. 
2001) implies binary nature. The ASCA X-ray spectra of 
ULXs have been described as the emission from optically 
thick accretion disks (Makishima et al. 2000 and refer- 
ences therein), however, the inferred disk temperatures, in 
the range of ~ 1 — 2 keV, are too high for BH masses im- 
plied by their X-ray luminosity. Recent observations with 
Chandra and XMM-Newton suggest lower temperatures in 
the range of ~ 100 — 500 eV (see review by Miller & Col- 
bert 2003 and references therein). Observations of spec- 
tral transitions between low (hard) and high (soft) state 
from two ULXs in IC 342 (Kubota et al. 2001) further 
demonstrate their similarity with the Galactic BH bina- 
ries. These observational facts suggest that ULXs may 
be more massive or scaled up versions of X-ray binaries. 
Recently, Ebisawa et al. (2003) argued that the high tem- 
peratures inferred from the X-ray spectra of some ULXs 
is high even for an accreting Kerr BH, and they explained 
the X-ray spectra in terms of a slim disk ( model without 
requiring intermediate mass BHs. 

Holmberg II X-l (hereafter HoII X-l) is nearby 
(3.39 Mpc; Karachentsev et al. 2002) and one of the 
brightest ULX in the sky. Its X-ray luminosity (Lx ~ 
10 40 erg s _1 ) corresponds to the Eddington luminosity of 



an 8OM BH. It resides near the edge of a dwarf star- 
forming galaxy and is embedded in an H II region. Also a 
compact Hell emission-line cloud has been discovered at the 
position of HoII X-l, indicating that this ULX is actually 
photo-ionizing the surrounding gas (Pakull & Mirioni 2001; 
Lehmann et al. 2004). The ROSAT HRI and PSPC obser- 
vations of HoII X-l revealed a point-like, variable source on 
scales of days and years (Zezas et al. 1999). The Chandra 
ACIS-S observations also revealed point-like X-ray source and 
an indication for a weak extended component (Lehmann et 
al 2004). Miyaji, Lehmann, & Hasinger (2001) presented the 
joint ROSAT and ASCA spectrum, which is well described by 
a power law (Fx ~ 1-9) and a soft excess component mod- 
eled either as an MCD (kT ~ 170 eV) or a mekal plasma 
(kT ~ 300 eV). In this paper, we present puzzling X-ray spec- 
tral variability of this well known ULX Holmberg II X-l ob- 
served with XMM-Newton. 

2. OBSERVATION AND DATA REDUCTION 

XMM-Newton observed Ho II X-l three times on 2002 
April 10, 16, and September 18 for 12.6, 13.8, and 6.9 ks, 
respectively. The EPIC PN (Striider et al. 2001) and MOS 
(Turner et al. 200) cameras were operated in full frame mode 
using the thin filter. The raw events were processed and filtered 
using the most recent updated calibration database and analy- 
sis software (SAS v5.4.l) available in December 2003. Clean- 
ing of the flaring particle background resulted in 'good' expo- 
sure times of 4.66, 3.95, and 4.25 ks for the EPIC PN obser- 
vations of 10 April, 16 April, and 18 September 2002, respec- 
tively. Events in the bad pixels and those adjacent pixels were 
discarded. Only events with pattern — 4 (single and double) 
for the PN and 0-12 for the MOS were selected. 

3. ANALYSIS & RESULTS 

We extracted full band light curves for HoII X-l from the PN 
data using circular regions of radii 40" centered at the source 
position. We also extracted the corresponding background light 
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curves from source free regions with the same bin sizes and ex- 
posure requirement. The background corrected PN light curves 
of HoII X-l are shown in Figure 1 for the three XMM-Newton 
observations. Small amplitude ( & 20%) fluctuations are seen 
in each of the three observations. A constant count rate fit 
to each light curve resulted in x 2 /dof = 62.0/51, 56.5/46, and 
40.3/46 for the observations of 10, 16 April 2002, and 18 Sep. 
2002. This suggests that X-ray emission was variable during 
the first two observations of April 2002. The mean flux level 
remained similar in these two observations. However, X-ray in- 
tensity from HoII X-l dropped by a factor of ~ 3 — 4 during the 
third observations of September 2002. For comparison we show 
the long-term light curve of HoII X-l in Fig. 1, derived for the 
common energy band of ROSAT, ASCA, and XMM-Newton. 
The ROSAT and ASCA fluxes were taken from Miyaji et al. 
(2001). It is clear that the source was in its brightest and 
faintest flux levels during the XMM-Newton observations of 10 
April 2002 and 18 Sep. 2002, respectively. 

Source spectra and the associated background spectra were 
extracted using similar extraction regions described above. The 
source spectra were grouped to a minimum of 20 counts and 
15 counts per channel for the 2002 April and September obser- 
vations, respectively. The spectra were analyzed using XSPEC 
11.3. All the errors quoted below were calculated at 90% con- 
fidence level for one interesting parameter. We found gener- 
ally good agreement between MOS and PN cameras in the 
0.3 — 10 keV bands. Therefore we present the spectral results 
obtained by fitting the same model jointly to the PN and MOS 
data while leaving the relative normalizations to vary. We fitted 
an absorbed power- law model to the spectra of 10 April 2002, 
this resulted in \ 2 = 933.4 for 758 degrees of freedom (dof). 
The best-fit absorption column is A^h ~ 1.9 x 10 21 cm~ 2 , much 
higher than the Galactic column of 3.4 x 10 20 cm~ 2 (Dickey 
& Lockman 1990) along the direction of Ho II X-l. The ex- 
cess column, inferred from an additional absorption model, is 
1.61JJ; 2 x 10 21 cm" 2 . Addition of blackbody (BB) model im- 
proves the fit significantly (model 1; A\ 2 = —62.5 for two 
additional parameters). The BB temperature is kT ~ 140 eV 
and power-law photon index is Fx ~ 2.6. The EPIC spectral 
data of 10 April 2002 and the best-fit model is shown in Fig- 
ure 2. Replacing the simple BB component by an MCD (model 
2), appropriate for a standard thin accretion disk, resulted in 
a slightly higher temperature. The best-fit model parameters 
are listed in Table 1. We followed the above steps to model 
the EPIC spectra of HoII X-l obtained on 16 April 2002 and 
18 Sep. 2002 and found that these spectra too are well de- 
scribed by the absorbed power law and BB model with the 
addition of the later component improving the fit significantly 
(Ax 2 = —61.8, and —26.7 for the spectra of April 10, and 
Sep. 18, respectively). The best-fit parameters for the three 
spectra are listed in Table 1. The X-ray spectrum of 18 Sep. 
2002 is weakest and steepest among the three spectra. The 
power-law normalization decreased by a factor of about four 
in five months. However, X-ray spectrum of HoII X-l is not 
the flattest at its brightest phase i.e., there is no evidence for 
correlation between the spectral slope and the flux. 

To illustrate the flux and spectral variability of HoII X-l, we 
compared each of the observed spectra with the best-fit spec- 
tral model to the spectrum of 10 April 2002. Figure 3 shows 
ratios of the individual spectra and the best-fit model spectrum 
of 10 April 2002. It is clear that the spectrum of HoII X-l did 
not change appreciably during the two observations in April 
2002, but became drastically weaker and steeper in Sep. 2002. 

The soft X-ray line- like features in the low state spectrum are 
suggestive of an optically-thin emission. Unfortunately, RGS 
data is too poor to constrain any line emission. The line-like 
feature at 2.5 keV, seen in the low state spectrum only (see 
Fig. 3), is detected at 2.8er level or 98% confidence level based 



on a maximum likelihood ratio (MLR) test (Ax 2 = -7.8 for 
two additional parameters) using the model 1 (see Table 1) as 
the continuum and a Gaussian line. This feature, also seen 
in the ASCA spectra (A X 2 = -7.1 for GIS+SIS data using 
the same continuum model), is likely the S XV Ka emission 
line. There may be some contribution from the radiative re- 
combination continuum (RRC) due to Si XIII. To investigate 
the effect of the presence of thin thermal emission on the spec- 
tral variability, we created a multi-component model consist- 
ing of a mekal plasma, a BB, and a power-law modified by 
the Galactic as well excess absorption. We fitted this model 
jointly to the PN spectra corresponding to the brightest flux 
state of 10 April 2002 and faintest flux state of 18 September 
2002. The parameters of the BB and the power-law compo- 
nents for the two data sets were varied independently, while 
the parameters of the mekal component for low state spectrum 
were tied to the corresponding parameters for the high state 
spectrum. This model resulted in a good fit (x 2 = 504.2 for 
499 dof). The unfolded spectra along with the models for the 
low and high states are plotted in Figure 4. The best-fit param- 
eters are: kT me k a i = 0.90^'g' 22 keV, emission measure EM = 
9.3 x 10 60 cm" 3 , kT BB {high) = U0+ 7 W cV, T(high) = 2.6±g;l, 
kT BB {low) = 106^3 eV, F(low) = 2.7t° 2 1 - The power-law 
normalization in the low state is a factor of 5.2 lower than 
that in the high state (2.4 x 10~ 3 photons cm" - 2 s- 1 keV" 1 at 
1 keV). The two photon indices are similar within errors. Thus, 
the softer spectrum in the low state can be explained as the 
drop of the flux in the power-law component without signifi- 
cantly changing its slope and a slightly softer BB component. 
The exclusion of the mekal component worsened the fit, result- 
ing in x = 512.6 for 501 dof. Thus, the presence of the mekal 
component is significant at a confidence level of 98.5% based 
on the MLR test. In the low state alone, the presence of the 
mekal component is significant at 99.3% level (Ax 2 = —10.1 
for two additional parameters). 

4. DISCUSSION 

The three XMM-Newton observations of 2002 have caught 
HoII X-l in its historical brightest and faintest flux levels 
(3.7 x 10~ 12 and 9.4 x 10" 13 erg cm" 2 s" 1 in the 0.5 - 2 keV 
band, respectively) ever observed with an X-ray satellite (see 
Fig. 1). Our spectral analysis of the three XMM-Newton obser- 
vations of HoII X-l confirm earlier R OS AT and ASCA results 
that the X-ray spectrum of HoII X-l consists of three spectral 
components - intrinsic absorption, a soft excess, and a power- 
law component. The intrinsic absorption column densities and 
temperatures (150 — 190 eV) of the soft component inferred 
from the MCD fits to three observations are similar to that 
obtained from the joint ROSAT and ASCA spectral analysis 
(Miyaji et al. 2001). However, the power-law photon indices 
are steeper (AF ~ 0.5 — 1.0) during the XMM-Newton obser- 
vations than that during the ASCA observations. 

A soft X-ray excess component, above a power-law contin- 
uum, is statistically required in all three observations and is 
well described by a simple BB or an MCD model. The best-fit 
temperatures are in the range 100 — 180 eV. Similar cool ther- 
mal components have been recently observed by Chandra and 
XMM-Newton from a number of bright ULXs e.g., M81 X-9 
(Miller, Fabian, & Miller 2003a), NGC 4038/4039 X-37 (Miller 
et al. 2003b), and NGC 5048 X-l (Karret et al. 2003). If the 
soft X-ray excess emission is attributed to the disk emission, as 
thought to be the case in BH X-ray binaries, the accretion disk 
of HoII X-l must be cool. In this scenario, the BB temperature 
provides a BH mass 

M ~ 10(fcT/1.2 keVy 4 {M/M Edd )M Q (1) 

(see e.g. Makishima et al. 2000). If we assume that a fraction 
e of the bolometric luminosity is emitted in the X-ray band 
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of 0.3 — 10 keV, then using the expression for the Eddington 
luminosity 



Lboi 



M 



L Edd M Edd £ x 1.28 x 10»( M/M ) 



(2) 



where Lx is in erg s 1 . Use of Eq. 1 with a BB temperature 
of kT = 170 eV for HoII X-l gives 



M 



Edd 



L x x 3.2 x 10- 42 



(e/0.1) 



(3) 



For L x ~ 2 x 10 40 erg s _1 as observed on 10 April 2002, 
M I MEdd ~ 0.2. If we adopt a BB temperature, kT = 110 eV, 
as suggested from the multi-component model including mekal 
plasma, the relative accretion rate is ~ 0.1. 

HoII X-l resides within a starburst region, therefore it is pos- 
sible that part of the X-ray emission arises from the starburst 
region or a supernova remnant. Our spectral modeling, after 
taking into account such a contribution as an optically-thin 
emission, suggests that the observed unusual spectral transi- 
tion to a low/soft state is due to the drop off the power-law flux 
alone without any change in the slope. A higher relative contri- 
bution of the optically-thin emission in the low state resulted 
in an overall softer spectrum. Thus, we suggest that the X-ray 
emission from HoII X-l consists of two distinct components: (i) 
a soft X-ray excess, described by an MCD (kT ~ 140 eV) and a 
power-law component (T ~ 2.7) similar to that observed from 



accreting BHs (AGNs and BH binaries) , and (it) a thin thermal 
component (kT ~ 0.9 keV, L(0.05-10 keV = 3xl0 38 erg s _1 ). 
In this case, the spectral variability of the accreting BH in HoII 
X-l is not unusual and similar to many BH X-ray binaries. 
Chandra ACIS-S observation barely resolved X-ray emission 
from HoII X-l, suggesting that the thermal component arises 
from a compact region (~ 10 pc; Lehmann et al. 2004). The es- 
timated luminosity of the thin thermal component arising from 
a compact region is not unusually large, as the cooling time, 
Tcooi ~ 50000 yrs is sufficiently long, where r coo ; was estimated 
from the temperature, electron density, size, and luminosity of 
the mekal plasma. A compact size and a high luminosity for the 
thin thermal component are suggestive of a young supernova 
remnant or hypernova remnant. The thin thermal emission ex- 
plains the line feature at ~ 2.5 keV, detected at 2.8a level in 
the low state spectrum only, as the S XV Ka line (see Figs. 3 
& 4), however, there may be some contribution from Si XIII 
RRC. High resolution X-ray observations will be required to 
investigate the nature of the soft component modeled here as 
thin thermal emission. 
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Fig. 1.— Left: EPIC PN light curves of Holmberg II X-l obtained from the three observations (a) 10 April 2002, (b) 16 April 2002, and 
(c) 18 September 2002. The light curves have been corrected for the background contribution and are shown with time bins of 100 s. Right: 
Long-term light curve of HoII X-l. ROSAT and ASCA fluxes were taken from Miyaji ct al. (2001). 
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Fig. 2. — The observed spectra of Holmberg II X-l and the best-fit model consisting of a power law and a BB modified by the Galactic as 
well as the intrinsic absorption. The bottom panel shows the deviations of the observed data points from the best-fit model in units of a. 
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Fig. 3.— A comparison of the observed EPIC PN data of (a) 10 April 2002, (b) 16 April 2002, and (c) 18 September 2002, and the 
best-fitting model of 10 April 2002 consisting of a power law and a BB modified by absorption. The ratios are shown for the normalization 
of the PN spectrum. The drop in flux below 0.5 keV in the spectrum of 10 April 2002 is present only in the PN data and is not seen in any 
of the MOS spectrum. 
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Fig. 4.— A comparison of the unfolded EPIC PN spectra of HoII X-l in the high (10 April 2002) and low (18 September 2002) flux states. 
In both cases, the best-fitting model consists of a BB, a power-law, and a mekal. The dotted lines are the model components (i) BB (low 
state), (ii) BB (high state), (Hi) mekal, (iv) power law. The data points with the solid line are the observed spectra and the best-fit model 
for (u) low state and (vi) high state. 
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Table 1 

Best- fit spectral model parameters for Holmberg II X-l. 



10 April 2002 16 April 2002 18 September 2002 

Parameter model l a model 2 a model l a model 2 a model l a model 2 a 



N H (10 21 cm" 2 ) 


i-6±8:? 


1.6^ 


1 4+0-3 

L -*-0.3 


1 r+0.3 
i '°-0.3 


1 4 +0.5 


1 5+°- 7 

1 -°-0.3 


fcT(eV) 


1411" 


187± 2 * 


128t 22 ; 




120t 22 . 


149t 28 7 


fbb b 


1.9 


2.7 


2.4 


3.4 




1.6 


T x 


2 fi4+ 06 

z -°^-0.06 


2 62+ - 07 

z - DZ -0.07 


2 40+ 07 

z -^ u -0.08 


9 07+0. 08 
z -° ' -0.07 


2.89+V.l? 


2.88±°;i 8 7 


fpL° 


10.1 


9.7 


6.9 


6.7 


2.4 


2.3 


f d 

J int 


14.8 


15.2 


12.5 


12.9 


3.9 


4.4 


T c 

i-'in 


2.0 


2.1 


1.7 


1.8 


0.5 


0.6 




870.7/756 


874.5/756 


606.9/562 


608.9/562 


270.0/268 


271.7/268 



"model 1 - simple BB and power-law model modified by the Galactic as well intrinsic absorption; model 2 - same as model 1 with the simple 
BB replaced with MCD. 
^Unabsorbed BB flux in units of 10 -12 erg cm -2 s^ 1 and in the 0.3 — 2 keV band. 
c Unabsorbed power-law flux in units of 10 — 12 erg cm -2 s -1 and in the 0.3 — 2 keV band. 
d Unabsorbcd flux in the 0.3 — 10 keV band and in the units of 10 12 erg cm -2 s — 1 . 
e 0.3 — 10 keV intrinsic luminosity in units of 10 40 erg s _1 assuming a distance of 3.39 Mpc. 



